Graphene via large N I: Renormalization 
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We analyze the competing effects of moderate to strong Coulomb electron-electron interactions 
and weak quenched disorder in graphene. Using a one-loop renormalization group calculation con- 
trolled within the large- iV approximation, we demonstrate that, at successively lower energy (tem- 
perature or chemical potential) scales, a type of non-Abelian vector potential disorder always asserts 
itself as the dominant elastic scattering mechanism for generic short-ranged microscopic defect dis- 
tributions. Vector potential disorder is tied to both elastic lattice deformations ("ripples") and 
topological lattice defects. We identify several well-defined scaling regimes, for which we provide 
scaling predictions for the electrical conductivity and thermopower, valid when the inelastic life- 
time due to interactions exceeds the elastic lifetime due to disorder. Coulomb interaction effects 
should figure strongly into the physics of suspended graphene films, where r s > 1; we expect vector 
potential disorder to play an important role in the description of transport in such films. 
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I. INTRODUCTION 

The recent experimental realization of graphene, 1-3 a 
single atomic monolayer of graphite, has refocused atten- 
tion upon many fundamental questions regarding elec- 
tronic transport in low dimensions. Its remarkable band- 
structure at zero doping has made graphene a candidate 
solid state analog for high energy particle phenomena, 
e.g. that described by the theories of quantum electrody- 
namics (QED) or chromodynamics (QCD). 

Given the relative simplicity of graphene as solid state 
system, it seemed not unreasonable to hope that a com- 
plete understanding of its basic electronic properties 
would emerge quickly. This expectation has been under- 
mined by several puzzling results in the first graphene 
experiments. 1-3 Perhaps most paradoxical is the "quasi- 
ballistic" nature of its electronic transport: at zero mag- 
netic field, graphene's conductivity is a linear, particle- 
hole symmetric function of carrier density that varies 
little over a temperature range of several orders of 
magnitude, 1-3 down to the lowest temperatures so far 
measured (30 mK). 4 At exactly zero doping, a condition 
termed the "Dirac point," the conductivity assumes a 
minimum value of order the conductance quantum e 2 /h; 
this "minimum metallic" conductivity varies only weakly 
between different samples possessing mobilities spanning 
an order of magnitude. 1-4 

One expects that any disorder present must be play- 
ing a crucial role in limiting the low-temperature con- 
ductivity of undoped graphene. Potentially important 
sources of disorder in the experiments performed in 
Rcfs. 1-5 include remote charged impurities in the Si02 
substrate, 6-10 as well as corrugations or "ripples" in 
the graphene sheet. 5,11-13 A semiclassical computation 
within the self-consistent Born approximation (SCBA) at 
the Dirac point gives a "bare" minimum metallic conduc- 
tivity independent of the disorder strength, and of order 



the conductance quantum, 



14-17 



m apparent agreement 



with experiment. Unfortunately, the SCBA is known to 



be inconsistent for massless Dirac fermions in 2D. 18,19 
Das Sarma ct al. 8,9 have employed an alternative mean 
field treatment focused upon charged impurity scatter- 
ing, and have concluded that the minimum metallic con- 
ductivity is in fact not universal, but should exhibit 
dependence upon the disorder distribution. Near zero 
doping, Cheianov et al. 10 have further argued that per- 
colation effects become important; these authors used 
semiclassical percolation theory to derive a scaling form 
of the conductivity in terms of the transmission coeffi- 
cient between electron and hole puddles. Other recent 
analytical 12, 19-22 and numerical 7,23,24 work has focused 
upon the localization physics of massless Dirac electrons 
in the presence of various types of disorder. These and 
previous 18, 25-29 studies have demonstrated that Dirac 
electrons may evade Anderson localization for certain 
fine-tuned disorder distributions, but the relevance of 
these results to graphene remains unclear. 

All of the references discussed above essentially ig- 
nore a potentially important aspect of graphene physics, 
that of electron-electron interactions. Carriers native to 
the carbon sheet arc ineffective at screening the long- 
ranged Coulomb potential for a system near zero dop- 
ing. The dimensionlcss strength of the Coulomb inter- 
actions r s r~j e 2 /ehvF, with vf the Fermi velocity and 
e the (effective) dielectric constant, is close to one in 
the experiments, 1-5 and should exceed two for suspended 
films. 11 This is not large by the standards of bulk 3D 
metals. However, it is important to stress that, near the 
Dirac point, Fermi liquid theory is not expected to apply, 
at least in the absence of disorder. Instead, the situation 
in graphene is analogous to that in the low-energy do- 
mains of QED and QCD in 3 + 1 spacctime dimensions; 
the smallness of the fine structure constant allows for 
a consistent perturbative expansion in QED for length 
scales of order or larger than the Compton wavelength. 
In the case of QCD, however, the coupling strength is 
not small at atomic distances, and in fact grows ever 
larger with increasing length scale. It is not apriori clear 
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that the Coulomb intcrparticlc interactions in graphene 
can be treated using perturbation theory in the coupling 
strength. 

For weak Coulomb interactions and disorder, one can 
use the perturbative renormalization group (RG) to de- 
termine whether disorder or interaction effects should be- 
come dominant at low temperatures. In this paper, we 
employ a large- AT generalization of the graphene field the- 
ory to treat the effects of moderate to strong Coulomb 
interactions; in physical graphene, N = 4 (see Sec. II, 
below). We study the resulting theory using a one- loop 
RG calculation perturbative in 1/N and in the disorder 
strength. 

The RG has been used many times before to study 
disorder effects in 2D Dirac fermions, 19,25 ' 27,28 ' 30 as well 
as the competition between disorder and interaction 
effects in such systems. 31-34 Using a non-interacting 
model of graphene, Aleiner and Efetov 19 showed that 
a particular type of disorder, associated with long- 
wavelength potential fluctuations consistent with, e.g., 
screened remote charged substrate impurities, dominates 
the flow to strong coupling in the absence of interactions. 
Stauber, Guinea, and Vozmcdiano 33 used the RG to treat 
both weak disorder and weak Coulomb interactions in 
graphene; they recovered results previously obtained by 
Ye and Sachdev. 31 ' 32 In particular, the weak-coupling RG 
for graphene demonstrates that a) in the absence of dis- 
order, Coulomb interactions are irrelevant (in the sense of 
the RG); the non- interacting, clean Dirac description is 
stable, and b) in the presence of generic disorder, the non- 
interacting, clean Dirac fixed point becomes unstable, and 
the system flows toward strong disorder and interaction 
coupling. For weak Coulomb interactions, the flow to 
strong coupling is again dominated by scalar potential 
disorder, consistent with the formation of local electron 
and hole "puddles." This picture would seem to sug- 
gest that the asymptotic low-energy physics in disordered 
graphene should be the same as that in an "ordinary," 
diffusive 2D electron gas (but see Refs. 21-24). Multiple 
crossover regimes are possible depending upon assump- 
tions regarding the microscopic disorder distribution. 19 ' 20 

Our large- N RG treatment of graphene follows closely 
a previous calculation of Ye, 31,32 who employed the same 
methodology to study strong Coulomb interaction and 
weak disorder effects in a Dirac system possessing a sin- 
gle Fermi point (generalized to N flavors). Different 
from Ye, we treat the case of two inequivalent Fermi 
points (valleys); each valley is generalized to N/2 identi- 
cal flavors. Our calculational framework therefore main- 
tains the important distinction between disorder poten- 
tials that scatter particles between inequivalent valleys 
and those that do not. Further, we correct a calcula- 
tional mistake made in Ref. 32, the correct determina- 
tion of which crucially affects our results. (See Sec. Ill A 
and Appendix A for details.) Son 35 recently revisited 
the problem of Coulomb interactions at large- N for clean 
Dirac electrons in graphene; in the limit of vanishing dis- 
order, we recover the findings of Ref. 35. 



We now summarize our primary results. We show that, 
within the large N approximation, disordered graphene 
resembles QCD more than QED. The theory becomes 
more strongly coupled at lower temperatures and longer 
length scales. This eventually leads to a breakdown of 
our RG calculation, so that we cannot determine the 
ground state phase within the framework employed here. 
However, we are able to identify various scaling regimes 
which may play a role in future experimental observa- 
tions, particularly in suspended graphene films. 11 Our 
main results are obtained via numerical integration of 
the RG flow equations, and are as follows: a) The non- 
interacting, non-disordered Dirac fixed point is unstable 
upon the incorporation of generic, arbitrarily weak dis- 
order, b) For the case of weak (r s C 1) or vanishing 
Coulomb interactions, the RG flows to strong scalar po- 
tential disorder, consistent with the electron-hole puddle 
picture. 8 ~ 10 > 19 c) Most importantly, for moderate (r s ~ 1) 
to strong (r s 3> 1) Coulomb interactions, scalar potential 
disorder fluctuations are parametrically cut off via screen- 
ing. The system flows to strong interaction and disorder 
coupling, but now, the runaway flow is dominated by 
different type of disorder, corresponding to a quenched 
SU(2) non-Abelian vector potential in the effective Dirac 
electron theory. This vector potential is implicated in 
ripples, 5,11 ' 12,20 as well as the representation of topolog- 
ical lattice defects, 36 in the low-energy theory. Finally, 
we use the RG to predict logarithmic temperature scal- 
ing in the dc conductivity within the Drudc/Boltzmann 
transport regime, and we also discuss thermal transport. 

The importance of quenched vector potential disor- 
der in the context of graphene has been stressed sev- 
eral times before in the literature. 5,11 ' 12,2 ° All of these 
previous works neglected interparticle interactions, how- 
ever, and overlooked the important fact that there is no 
physically plausible way of disordering a graphene sheet 
that would produce only non-Abelian vector potential 
disorder in isolation (see Sec. II A, below). Our result 
establishes that, for moderate to strong Coulomb inter- 
actions in the large-A^ approximation, non-Abelian vec- 
tor potential disorder emerges generically as the dom- 
inant source of randomness. Finally, it is intriguing 
to note that the non-interacting graphene theory sub- 
ject only to non-Abelian vector potential disorder con- 
stitutes a very unusual, manifestly particle-hole sym- 
metric quantum disorder model, for which several ex- 
act and/or non-perturbative results are known. 18 ' 20 ' 26 
This non-Abelian Dirac disorder model is predicted to 
avoid Anderson localization, possessing instead a pertur- 
batively inaccessible, critical ground state, 26 a vanishing 
disorder-averaged, single-particle density of states, 18 and 
a non-zero conductivity of order e 2 /h, independent of the 
disorder strength. 20 

The outline of this paper is as follows. In Sec. II, we 
write down the low energy field theory for disordered, in- 
teracting graphene, and generalize it to the case of largc- 
N in order to treat strong Coulomb interactions. We dis- 
cuss the microscopic interpretation of the various types of 
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disorder that appear in the low-energy theory, and we re- 
view the quantum disorder universality classes obtained 
by fine-tuning this distribution. In Sec. Ill, we perform 
a one-loop renormalization group calculation upon the 
model defined in Sec. II. In Sec. IV, we present the results 
of this calculation, which consist of flow equations for 
the various parameters defining the graphene field the- 
ory. We analyze these equations and attempt to gleam 
information about graphene's phase diagram in an ab- 
stract disorder-interaction coupling constant space. Ex- 
perimental consequences of our results in the form of scal- 
ing predictions for the dc conductivity and thermopower 
are summarized and discussed in Sec. V. 

A variety of elaborations and extensions are relegated 
to the Appendices. Some technical details of the calcu- 
lation presented in Sec. Ill are relegated to Appendix 
A, while a survey of several (unstable) fixed line struc- 
tures appearing in the flow equations stated in Sec. IV is 
presented in Appendix B. The scaling predictions enu- 
merated in Sec. V arc derived in Appendix C. 



II. MODEL AND SETUP 

The plan of this section is as follows: In subsection 
II A, we write down an effective field theory for graphene 
incorporating both quenched disorder and Coulomb in- 
terparticle interactions, and we discuss (1) the role of dis- 
crete symmetries and (2) the different universality classes 
of disordered quantum systems that arise by fine-tuning 
the details of the impurity or defect distribution. Most 
of this material is not new, see e.g. Refs. 19,20; the pur- 
pose of this section is to provide a context for the re- 
sults discussed in Sec. IV, and to establish notation. In 
subsection II B, we construct a large- N generalization of 
our graphene model, suitable for studying the effects of 
strong Coulomb interactions. 



A. Lattice model and effective field theory 

The bandstructure of undoped graphene is well- 
approximated by the tight-binding Hamiltonian 



Hn = -t 



E c L(r)c Bs (r')+H.c. 



(2.1) 



(rr'}, 



where c^ s (r) and c b - s {t') are creation and annihilation 
operators for electrons on the A and B sublattices of 
the bipartite honeycomb lattice, respectively. The index 
s G {Til} denotes spin-1/2 components of the lattice 
electrons. The hopping amplitude t in Eq. (2.1) is purely 
real, and the sum on (rr') runs over all nearest-neighbor 
A — B bonds in the graphene sheet. (Further- neighbor 
hopping may be incorporated as a perturbation). Hence- 
forth we adopt units such that H = ks = 1- 

A low-energy effective field theory for undoped 
graphene obtains by linearizing the bandstructure of 



Eq. (2.1) in the vicinity of the two inequivalcnt K and K' 
Fermi points (also termed "valleys"), where the particle- 
hole symmetric energy bands meet. Retaining only low- 
energy modes, we write 



5 (r) 



- k - r V£'(r), 



(2.2) 



where ±kp locates the K and K' points, respectively, 
and a G {^4, B} labels the sublattice species. 
We assemble the Dirac spinor 



^(k) = 



<(k)' 
<(k) 



(2.3) 



(The spin index s has been suppressed in this equation.) 
In the low-energy theory, we will consider states annihi- 
lated by V'(k) with momentum k < A <C kp, where A is 
a hard cutoff. 

The Dirac spinor defined by Eq. (2.3) is an eight- 
component object, V - * ipa: with index a G {1, . . . , 8}. 
The eight components arise from the direct product 
of indices in the 2-dimensional sublattice {A,B}, val- 
ley {K, K'}, and spin 1/2 {|, 1} component subspaces. 
Later, we will generalize ip a to 2N components in order 
to perform an expansion in 1/N. We introduce two com- 
muting sets of Pauli matrices: the matrix a a acts in the 
sublattice {A, B} space of Eq. (2.3), while the matrix hp 
acts in the valley {K, K'} space, with a,/3& {1, 2, 3}. 37 

Below, we write down the continuum field theory for 
graphene that incorporates the effects of both quenched 
disorder and Coulomb electron-electron interactions [see 
Eq. (2.8)]. The clean model in Eq. (2.1) possesses three 
crucial symmetries that might plausibly survive the in- 
corporation of disorder: particle-hole symmetry (PH), 
time-reversal invariance (TRI), and spin SU(2) rotational 
symmetry. In the low-energy theory, the former are en- 
coded in the operator level transformations 

V>(r)->-* 1 K 1 [V>(r)t] T , (PH); (2.4) 

Hr) ^a 2 k 2 ^(r), (TRI). (2.5) 

In this equation, T denotes the ordinary (matrix) trans- 
pose. The particle-hole (time-reversal) transformation 
defined by these equations is unitary (antiunitary). 
[Schematically writing H = ip^htp, c.f. Eq. (2.8), Eq. (2.5) 
implies that TRI imposes the condition & 2 K2 h* a 2 £2 = h 
upon the single-particle Hamiltonian h.] 3S The product 
of the operations given by Eqs. (2.4) and (2.5) is a so- 
called "chiral" transformation, 28,39-42 and was denoted 
by the symbol "C z " in Ref. 20; in this same reference, 
the physical TRI [Eq. (2.5)] operation was labeled "To-" 
In addition, it is useful to define two additional trans- 
formations, which do not correspond to microscopic sym- 
metry operations in the tight-binding honeycomb model 
of graphene, but nevertheless may play important roles as 
emergent symmetries at a special critical point, or when 
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additional restrictions arc placed upon the disorder dis- 
tribution: 



lfr(r) 



-a 1 ^ [^(r) t 
a 2 ^(r), 



(PH*); 
(TRI*). 



(2.6) 
(2.7) 



zero, while we may assign up to five independent param- 
eters {g u , gA, gA3, grn, 9v} to characterize their statistical 
fluctuations: 19 



The transformation PH* is defined to be unitary, while 
TRI* is antiunitary. The product of PH* [Eqs. (2.6)] and 
the physical TRI operation [Eq. (2.5)] is another type of 
"crural" transformation, denoted by the symbol "Co" in 
Rcf. 20; the TRI* operation [Eq. (2.7)] was labeled "T x " 
in this reference. 

The effective field theory for undoped graphene is en- 
capsulated by the Hamiltonian 

H= J d 2 r^{r) -iv F a ■ V + V(r) ^(r) 

+ w r Mrf fm^. (2 . 8) 

2 J \r- r'| 

In this equation, vp ~ 3t/2 is the Fermi velocity, V(r) 
is a matrix-valued single-particle potential that encodes 
the structure of the quenched disorder (discussed below), 
and W ~ e 2 /e is the microscopic strength of the Coulomb 
interactions. 

We assume that both TRI [Eq. (2.5)] and spin SU(2) 
rotational symmetry survive in every static realization 
of disorder. Our study thus excludes the effects of ex- 
ternal magnetic fields or spin-flip impurities, but is ap- 
propriate to remote charged impurities, "ripples," 5,11-13 
non-magnetic interstitial atoms, and topological defects 
such as dislocations and bound disclination pairs. Wc 
will be interested in the properties of a graphene system 
ensemble-averaged over many realizations of disorder; we 
presume that such averaging restores invariance under 
PH transformations [Eq. (2.4)], as well as all honeycomb 
lattice space group symmetry operations (translations, 
rotations, reflections). Then the disorder potential V(r) 
appearing in Eq. (2.8) is conveniently parameterized as 

V(r) =u(r) i + A|(r) <r% + A%(t) & & k 3 

+ mP(r)a 3 kp + v(r)a 3 k 3 . (2.9) 

In this equation, the barred indices a and f3 are under- 
stood to be summed over the "spatial" Pauli matrix com- 
ponents 



Me {1,2}. 



(2.10) 



Note that Eq. (2.9) implies that potentials {A^} and 
{m* 9 } scatter between inequivalent Fermi nodes (i.e. cou- 
ple to the valley space Pauli matrices h\ or fa®); u, 
{A?}, and v do not. Imposing statistical invariance un- 
der particle-hole [Eq. (2.4)] and honeycomb lattice space 
group transformations, and assuming white noise (short- 
ranged), Gaussian-correlated disorder, we must set the 
average value of all disorder potentials in Eq. (2.9) to 



u(r) u(r') 


= 2ng u v 2 F S^(r-r'), 


(2.11a) 


4(04(r') 


= 2irg A v 2 F 6 &rx 6^6W(r-r'), 


(2.11b) 


A|(r)4(r') 


= 27rg A3 v 2 F 6 arx S^(r-v'), 


(2.11c) 


777,0 (r) m^(r') 


= 2-Kg m v 2 F ^5^{v-v'), 


(2. lid) 


v(r) v(r') 


= 2ng v v 2 F S^(r-r'). 


(2.11c) 



Let us now take a brief excursion from setting up our 
model to review the physical interpretation of these var- 
ious disorder types described by Eqs. (2.9) and (2.11); 
see also Refs. 19,20,36. Wc will emphasize the connec- 
tion between the particular types of disorder, the "sym- 
metry" operations defined by Eqs. (2.4)-(2.7), and the 
different universality classes 39-41 of disordered quantum 
systems realized by fine-tuning the disorder distribution 
in graphene. 



1. Scalar potential disorder and the symplectic class All 

Remote charged impurities at or near the substrate 
surface may provide the dominant scattering mechanism 
for electrons in the experiments of Refs. 1-5. One ex- 
pects the potential arising from these impurities to vary 
relatively slowly on the graphene lattice scale, and there- 
fore to manifest itself primarily through the effective 
scalar disorder potential u(r) in the low-energy theory 
[Eqs. (2.8)— (2.11)]. We emphasize that the translation 
of any lattice-scale disorder potential into the single ef- 
fective disorder field it(r) is typically impossible, and 
all disorder types consistent with the underlying sym- 
metry [i.e. time-reversal invariance, Eq. (2.5)] will ap- 
pear in the low-energy theory. Nevertheless, theoreti- 
cally we are free to fine-tune the disorder distribution in 
Eqs. (2.8), (2.9), and (2.11) by setting to zero all disor- 
der parameters except g u , i.e. only u(r) nonzero. The re- 
sulting theory possesses a large SU(4) symmetry, present 
in every fixed disorder realization, associated with uni- 
tary transformations in (valley) <S> (spin-1/2) space. In 
every realization of u(r), this theory is invariant under 
the alternative "time-reversal" transformation TRI* , de- 
fined by Eq. (2.7). The advent of TRI* places the non- 
interacting system [Eqs. (2.8) and (2.9) with W = and 
only u(r) ^ 0] into the "symplectic" or spin-orbit uni- 
versality class All of disordered metals, 19-25 rather than 
the orthogonal class AI nominally expected on the ba- 
sis of spin-1/2 SU(2) rotational symmetry and "physical" 
TRI [Eq. (2. 5)]. 43 (We have adopted the nomenclature for 
quantum disorder classes employed in Refs. 39-41.) Un- 
like the orthogonal class, the symplectic class is known to 
possess (in the absence of interactions) a disorder-driven 
metal-insulator transition in 2D. 44 The (de-)localization 
of Dirac electrons subject only to scalar potential disor- 
der has been studied in Refs. 21-24. 
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2. Vector potential disorder and ripples, topological defects, 
CI and BDI classes 

In the context of graphene, the intravallcy Abclian 
vector potential {A^} in Eqs. (2.8) and (2.9) arises 
in the description of elastic lattice deformations or 
"ripples." 5,11 ' 12 ' 2 ° Ripples constitute a potentially im- 
portant source of scattering in the electronic transport 
of graphene, the presence of which has been directly 
observed and indirectly inferred in substrate-supported 5 
and suspended 11 graphene samples, respectively. Topo- 
logical lattice defects constitute another class of sample 
imperfections. In graphene, lattice dislocations can be 
modeled by point flux insertions of the potential {^4^} in 
the low-energy Dirac theory [Eq. (2.8)], while the descrip- 
tion of lattice disclinations requires the intervalley com- 
ponents {A?}. 3e Taken together, {A^} and {A%} form 
a quenched, non-Abelian SU(2) vector potential realized 
in valley space. 

The graphene field theory [Eqs. (2.8) and (2.9)] pos- 
sessing only quenched SU(2) vector potential disorder, 
obtained by fine-tuning u(r) = mr{r) = v(r) = 0, 
is invariant (in every disorder realization) under the 
alternative "particle- hole" transformation PH*, defined 
by Eq. (2.6). In the absence of interparticle interac- 
tions [W = in Eq. (2.8)], the corresponding disorder- 
averaged Dirac model with g u = g rn = g v = in 
Eq. (2.11) is a critical (i.e. not Anderson localized) 
conformal field theory (CFT), for which several exact 
and/or non-perturbative results are known. 18 ' 26 This 
CFT realizes a (variant of) the non-Wigner-Dyson uni- 
versality class CI, typically associated with disordered 
superconductors. 18 ' 40,41 Intriguingly, the Dirac theory 
with only non-Abelian disorder is known to possess a 
vanishing density of states, 18 and is predicted to ex- 
hibit a nonzero conductance, independent of the disor- 
der strength. 20 We note, however, that the PH* sym- 
metry transformation given by Eq. (2.6) does not cor- 
respond to any local operation in terms of the lattice 
electrons appearing in the model defined by Eq. (2.1). 
As a consequence, neither elastic deformations nor topo- 
logical defects occuring in the lattice-scale description of 
graphene translate into pure non-Abelian gauge disorder 
in the low-energy Dirac theory. For example, the field 
theory description of lattice dislocations and disclinations 
incorporates the other disorder types {u(r), m^(r), v(r)} 
into the cores of the defects, so that an idealized graphene 
sheet subject only to topological disorder does not realize 
the CI class CFT. 

Potentials {A^} and {m^} preserve both TRI and 
(physical) particle-hole symmetry (PH) [Eq. (2.4)]. 
These particle-hole symmetric potentials can be mapped 
to a particular type of "microscopic" disorder consisting 
of real, nearest-neighbor random hopping on the honey- 
comb lattice. Such disorder might be realized in principle 
in a graphene sheet by introducing lattice vacancies, if 
the only effect of removing a carbon atom is to sever the 
3 nearest-neighbor electronic bonds surrounding the va- 



cancy (i.e., treating the vacancy in the unitary limit). In 
the presence of TRI [Eq. (2.5)] and spin SU(2) rotational 
symmetry, and the absence of interactions (w = 0) and 
further-neighbor hopping, such a model with only gA3 
and g rn non-zero resides in the "crural" symmetry class 
BDI, 27 ' 28,41 and is another example of a critical, delocal- 
ized disordered-Dirac model in 2D for which exact and/or 
non-perturbative results arc available. 28 



3. Staggered on-site potential disorder and class D 



Finally, the disorder potential v(r), which couples to 
a matrix diagonal in both sublattice (a) and valley (k) 
spaces [Eq. (2.9)], may arise in concert with u(r), e.g., in 
the description of interstitial impurities. The imposition 
of a uniform v(r) = vo potential is consistent with the 
application of an on-site, A-B sublatticc-staggcred chem- 
ical potential to the underlying honeycomb lattice model 
[Eq. (2.1)]. The non-interacting Dirac theory with only 
staggered potential disorder, Eqs. (2.8) and (2.11) with 
W = and only g v > 0, realizes the non-Wigner-Dyson 
class D, and is formally related to the random bond Ising 
model in 2D. 25 < 45 - 48 

We see that by fine-tuning the details of the disor- 
der distribution within the confines of a graphene sheet 
that possesses, on average, the full space group symme- 
try of the honeycomb lattice, several universality classes 
of (non-interacting) disordered quantum systems may be 
realized. The realization of these classes is essentially 
equivalent to imposing various combinations of the dis- 
crete symmetries defined by Eqs. (2.4)-(2.7). We em- 
phasize, however, that in the current situation of zero- 
magnetic field graphene experiments, 1-5 only the physical 
time-reversal symmetry [Eq. (2.5)] presumably persists, 
and the ultimate low-energy, asymptotic physics must be 
described by the orthogonal Wigner-Dyson class AI. 19 
Moreover, the Coulomb interactions are expected to play 
a strong role, as we will demonstrate in this paper. 

We now complete the formal setup of our graphene 
field theory. We employ the zero-temperature, imagi- 
nary time path integral formalism, implementing repli- 
cas ip — > ip l a in order to average over the disorder. (The 
index i G {1, 2, . . . , n}, with n — ► at the end of the cal- 
culation). The post-disorder averaged, imaginary time, 
coherent state path integral is 



Z= \ V^V^e~ S , (2.12) 
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with 

S = J drd 2 rft (8r - iv F & • V) ^ 



2nG»v F 



drdr' d 2 r^^(r) {t') 



dr dr' dfv^a^h^^r) ip a* 1 k^ 1 (t '), 

(2.13) 



where summation is implied over the repeated replica 
and disorder vertex {/i, v} indices (defined below). 
The disorder is encoded in the terms on the third and 
fourth lines of Eq. (2.13), the latter of which involves the 
disorder "metric" 



9a 9a 9A3 
9a 9a 9A3 

9m 9m 9v 



(2.14) 



[c.f. Eqs. (2.9) and (2.11)]. couples to the direct prod- 
uct o^Ky <g> o^ky, with n, v S {1, 2, 3}. 

Assuming that frequency carries the (generally scale- 
dependent) "engineering" dimension 

kl = ( 2 - 15 ) 

measured in inverse-length units, we assign the dimension 
of the %j) field as 



[V(r,r)] = l. 



(2.16) 



Then wc find that the Fermi velocity v F and the Coulomb 
interaction strength VV share the dimension 

[v F ] = [W] = z - 1, (2.17) 

while the disorder strengths {<?u,!?^} are dimensionless. 



B. 



The theory at TV = oo 



The field theory given by Eqs. (2.12)-(2.14) was stud- 
ied by Stauber et al. in Ref. 33 for the case of weak dis- 
order and weak Coulomb interactions. (Their results are 
discussed in Sec. IV.) In graphene, the effective Coulomb 
interaction strength is of order unity, calling into ques- 
tion the usefulness of a perturbative expansion in powers 
of W [Eq. (2.13)]. In this paper, we employ a large- 
TV generalization in order to treat the case of moderate 
to strong Coulomb interactions. The expansion param- 
eter 1/TV = 1/4 for physical graphene, which is at least 
smaller than the bare dimensionless Coulomb interaction 
strength. Our treatment follows closely that of Ye in 
Ref. 32, wherein the author considered a large- TV gener- 
alization of a Dirac theory possessing a single valley. 



D 11 ib 1 ia 



D i.2 "VWW' 




/ a 



ib 



1° 



\jd 



FIG. 1: Diagrammatic elements of the bare Feynman rules. 
The associated amplitudes are summarized in Table I. 



We generalize the Dirac spinor ip l a [Eq. (2.3)] to 2TV 3> 
1 components in the index a, with physical graphene cor- 
responding to TV = 4. [i is the replica index, defined 
above Eq. (2.12).] Crucially, we will consider only even 
TV, retaining the partioning of TV into two inequivalent 
Fermi point sectors; as above, we will continue to use 
the Pauli matrices k m , m G {1,2,3}, to address compo- 
nents in this valley space. In other words, we may write 
TV = 2(2,s + 1), where we take the spin s ^> 1 — formally, 
this is a large spin expansion. At large TV, we consider 
only the spin- independent, "physical" disorder potentials 
defined by Eqs. (2.9) and (2.11), acting in the product of 
sublattice (a) and valley (k) spaces. 

We begin by decoupling the Coulomb interaction in 
Eqs. (2.12) and (2.13) with a Hubbard-Stratonovich 
transformation, using the temporal gauge field a 1 (i is 
again the replica index): 



VipVipVa e 



-s 



where 



d T -i\/ — a l (r,0) - iv F ar • V 



S -> / drd 2 r^ 1 

i J drd 2 rdz [[Va i (r,z)] 2 + [d z a l (r, z)] 2 ] 



(2.18) 



mg u vp 



drdr d z v^^{r) ^^{t ') 

drdr 1 d 2 x^ l b»h v i> l {T)^G ll k v ^(T l ). 

(2.19) 



( 




FIG. 2: The fermion propagator at TV = oo. 
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( WW 1 - n a ) _1 



n 



FIG. 3: The gauge (Coulomb) propagator at N = 



The temporal gauge field a 1 resides in the full 3 + 1 space- 
time (t, r, z), while the grapheme field theory is confined 
to the plane 2 = 0. As a result, the coefficient in front of 
the kinetic term on the second line of Eq. (2.19) cannot 
be renormalized, as we check explicitly in Sec. IIIC. In 
Eq. (2.19), we have replaced W — > w/2irN; we keep w 
finite and non-zero so as to obtain a well-defined theory 
in the limit TV — > oo. The bare Feynman rules are stated 
in Fig. 1 and in Table I. In this paper, we denote each 
diagram and its corresponding amplitude by the same 
symbol, D,„.„, where m is the figure number and n is the 
particular diagram in question. 

The TV — > oo theory, in the absence of disorder, is 
completely characterized by the one- and two-particle 
self-energies depicted in Figs. 2 and 3, respectively. For 
the special case of massless Dirac fermions studied here, 
however, particle-hole symmetry prevents the generation 
of a chemical potential shift, so that the electronic self- 
energy Ily, vanishes at TV = oo (and zero disorder), and 
the large- TV fermion propagator is the same as Di.i (Fig. 1 
and Table I). 

The gauge field self-energy at TV = oo is depicted in 
Fig. 3; it is given by the expression 



n a (ft„,q) = (i^jN) 2 V {n n ,q), 



(2.20) 



where the Dirac fermion polarization bubble T>o(fl n ,q) 
is 

TV o 2 

Z>o(fi„,q) = --=====. (2.21) 



16 y/W+vj.q* 

Incorporating the self-energy from Eq. (2.20), the 
propagator for the temporal gauge field at TV = oo and 



TABLE I: Factors associated to bare propagators and vertices 
depicted in Fig. 1. 



Di.i 

Dl.2 
Dl.3 

Dl.4 



(iuJn + VF<T ■ lt)a,b 



U>n + VFk 2 



2%g u v F 5 a ,b5 c ,d 
+2Tvg^v F {a fJ 'k h ,) aib (a^k,) Ctd 



zero disorder is 



D 



3 — ~ 



q ^v'j-q 2 + ni + wv F \q\ 



(2.22) 



where we have introduced the dimcnsionless Coulomb in- 
teraction strength 



w 
16vf 



(2.23) 



[c.f. Eq. (2.17)]. For the physical case of TV = 4, the di- 
mensionless parameter w is essentially the Coulomb in- 
teraction constant r s : 



w 



7r e 
2 evF 



(2.24) 



The scalar potential disorder, characterized by the dis- 
order strength g u [Eqs. (2.9) and (2.11), Fig. 1 and Ta- 
ble I] is parametrically screened by the Dirac sea in the 
large TV limit; the screening affects only the scalar poten- 
tial disorder, since the other disorder fields in Eq. (2.11) 
can be understood as arising from charge-neutral de- 
fects that produce potentials lacking a monopole (i.e. 
lowest-order multipole) moment. The diagram D4 pic- 
tured in Fig. 4 gives the complete dressed disorder vertex 
at TV = 00: 



D 4 =2ng u v F 5 a . b S c . d 



(2.25) 



where we have defined the screened scalar potential dis- 
order strength 



9u 



flu 



(1 + w) 2 



(2.26) 



In the large- TV theory, it is the screened strength g~ u , 
rather than the bare parameter g u , that appears in phys- 
ical observables (Sec. V and Appendix C), and there- 
fore characterizes the effective strength of scalar potential 
fluctuations. 

In summary, the Feynman rules for the large- TV ver- 
sion of graphene are given by Di 1 and D1.3 (Fig. 1 and 
Table I), and by D 3 and D 4 [Figs. 3 and 4, Eqs. (2.22) 
and (2.25), respectively]. 

i a j c 



ib jd 



FIG. 4: Screening of the scalar potential disorder g u due to 
Coulomb interactions in the TV = 00 limit. 



8 



D 



51 a) I 

<D„+<D 

k + I 



D 



5.2 



l-k 



FIG. 5: Corrections to the electronic self-energy at 



III. PERTURB ATIVE EXPANSION 

In this section, we perform a one-loop renormalization 
group (RG) calculation within the large- N generalization 
of graphene defined by Eqs. (2.18) and (2.19). Using the 
results derived below, we obtain flow equations for the 
disorder and interaction coupling strengths defined by 
Eqs. (2.8)-(2.11). These flow equations are stated and 
analyzed in Sec. IV. 



A. Electronic self-energy 

We begin with the renormalization of the electronic 
self-energy due to the disorder and Coulomb interactions 
(at finite, but large N). The diagrams are depicted in 
Fig. 5. We consider first the 1/N correction due to the 
interaction, diagram D5.1. Expanding in powers of the 
external frequency uj n and momentum k, we find that 



D 



5.1 



~r\ [wj x {w)iuj n +wf 2 (w)v F & ■ k]lnA, (3.1) 



where w is the dimcnsionlcss Coulomb interaction 
strength defined by Eq. (2.23), A is a hard momentum 
cutoff, and we have introduced 



V 



ttN' 



(3.2) 



In Eq. (3.1), the functions /i(u>) and fi{w) denote the 
following analytic functions of the interaction strength w 



fi(w) =- 



1 1 



i-vT 



2 VT 



w 2 \ 1 + Vl — w 2 



vT 



; arcsin(w) — w 



h(w) =- 



21 + Vi-w 2 



+ 



1 



7TW 



\/ 1 — w 2 arcsin(u>) — w 



(3.3a) 



(3.3b) 



It will prove useful to consider the sum 
h{w) =h{w) + f 2 (w) 



1 



2w 2 I 7T 



2 r arccos(iif) 

w + , ' - 1 



VT 



(3.4) 



The electronic self-energy contribution detailed in 
Eqs. (3.1)-(3.4) has been considered several times before 
in the literature, first by Ye and Sachdev, 31,32 and soon 
after by Gonzalez, Guinea, and Vozmediano. 49 The in- 
teraction self-energy at order 1 /N was revisited recently 
by Son. 35 Our results agree with those of Refs. 35,49, but 
disagree with those of Refs. 31,32. We derive Eqs. (3.1)- 
(3.4) explicitly in Appendix A. 

Next, we turn to the disorder diagram D5.2. We find 



D5.2 ~ Q s iu n In A, 



(3.5) 



where 



Q s =gu + % 



4g A + 2g A3 + 2g m + g v 



(3.6) 



[g u is the screened scalar potential disorder strength, de- 
fined by Eq. (2.26)]. 

The bare irreducible two-point vertex function is 



iT 



(0) 



Vp<7 



D 



5.1 



D 



5.27 



and the renormalization condition is 5C 



din A 



(V;r (0)N l 



0. 



(3.7) 



(3.8) 



where is the wavefunction renormalization of the ip 
field. Using Eqs. (3.1) and (3.5) in Eq. (3.8), and combin- 
ing the result with the dimensional analysis [Eqs. (2.15)- 
(2.17)], we obtain the flow equations 



dlnZ. 



■4, 



dl 
din vp 



dl 



-rjwfxiw) 

z-l-G s 



rjwf 3 (w) 



(3.9) 
(3.10) 



where I = — In A is the log of the RG length scale, and 
z is the (as yet unspecified) "dynamic critical exponent" 
[Eq. (2.15)]. 

Although we are employing a field-theoretic renormal- 
ization scheme, Eqs. (3.9) and (3.10) can be equivalcntly 
understood in a Wilsonian framework: after an integra- 
tion of high-momentum modes, one acquires corrections 
akin to those expressed in Eqs. (3.1) and (3.5). [Consult 
the action given by Eq. (2.13).] Rescaling momentum 
k — > k/b, frequency uj n — > uj n /b z , and field tp — > 1 ^ 2 ijj, 
one recovers Eqs. (3.9) and (3.10), where 6^1 



dl. 



B. Coulomb vertex 

One-loop corrections to the Coulomb vertex D1.3 
(Fig. 1 and Table I) are depicted in Fig. 6. Most of 
these diagrams involve a "3-electron ring," that is, a 
closed electron loop with three nodes. Such structures 
are generic to large- N expansions. We note that the 
graphene (masslcss Dirac) theory studied here possesses 
two simplifying features: 
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Y-'-V-S D 6.5 




FIG. 6: Coulomb vertex corrections at O (j?, {5u,{?^}) 



a) The sum of two counter propagating, but other- 
wise identical, m-electron rings vanishes for odd m 
whenever all m bilinear insertions involve the iden- 
tity matrix 1; see Fig. 7a. This result is simply a 
consequence of the odd parity nature of the mass- 
less Dirac Green's function G(r): 



G(r-r') = -G(r'-r). 



(3.11) 



b) An m-electron ring involving m — 1 identity matrix 
insertions and a single disorder matrix insertion of 
the form = a^k v vanishes trivially by the trace 
over valley (k) space; see Fig. 7b. 

As a result of this simplification, only diagrams Dg.i 
and D6.2 contribute. We find that 



D 



(j.i 



— ) rfwfi(w)bxA, 



and 



D 



6.2 



10 




AT 



G s In A. 



(3.12) 



(3.13) 



In these equations, fi(w) was defined by Eq. (3.3a), 
7] by Eq. (3.2), and Q s by Eq. (3.6). The evalua- 
tion of Eq. (3.12) closely parallels that of the electronic 
self-energy, the calculation of which is detailed in Ap- 
pendix A. 

The bare vertex function is 



iT 



(0) 

aijjijj 



1 + In A 



Qs-nwh^)]) ■ (3.14) 



Now, we need to know the wavefunction rcnormalizations 
for both the Dirac (ip) and gauge (a) fields to proceed; 



the latter requires the evaluation of the gauge field self- 
energy to one loop. Gauge invariance in fact requires 
that 



<iln Z a 
dm A 



0. 



We will check Eq. (3.15) in Sec. IIIC, below. 
The rcnormalization condition is therefore 



din Z, 



din. A din. A 



iT {0) -0 



(3.15) 



(3.16) 



which implies the following RG flow equation for the di- 
mcnsionless Coulomb interaction strength w [introduced 
in Eq. (2.23)]: 



dim 

~dT 



T]wfs(w) 



(3.17) 



[I = — In A in this equation]. To obtain Eq. (3.17), we 
have used Eqs. (3.4), (3.9), and (3.10). 



C. Gauge field self-energy 

The diagrammatic corrections to the gauge field self- 
energy to are depicted in Figs. 8 and 9. Because of the 
properties explicated in Fig. 7, the sum of the diagrams 
in Fig. 9 gives zero. 

All of the diagrams in Fig. 8 are formally two-loop, 
but nevertheless give only contributions to first order in 
In A. In other words, one of the two loops in each of these 
diagrams is always ultraviolet-finite, and we can obtain 
their contributions to the rcnormalization group using 
the hard momentum cutoff scheme implemented above. 

We begin with Ds.i and Ds.2, which give identical con- 
tributions. One finds that 



D 



8.2 



Vpq 2n2 +2Bw^v F j(n n ,<i), 



(3.18) 



a) 





1^3 =0 



b) 



= 



FIG. 7: a) Cancelation between counter-propagating 3- 
electron rings involving only identity matrix (1) insertions 
into the trace, b) 3-electron rings involving two identity ma- 
trix insertions and one disorder matrix of the form M pl , = 
a^kv vanish identically. 
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FIG. 8: Corrections to the gauge (Coulomb) field self-energy 
to 0(i {£1,^}). 



D 10.1 r a 






FIG. 10: Evaluation of "two-loop" diagrams Dg.s and Dg.e 
in Fig. 8. These graphs can be expressed in terms of prod- 
ucts of an ultraviolet-convergent polarization bubble and an 
ultraviolet-divergent vertex correction (see Fig. 6). 



where we have used Eq. (2.21), and the log-divergent 
constants A and B follow from Eq. (3.1): 



A = -^nw[h{w) - h(w)} In A, 



B = --r)wf 3 (w) In A. 



(3.19a) 
(3.19b) 



In Eq. (3.18), we have defined 

' [i{u + n n ) + v F a--(i + q)Y 

x [iuj + t;^<T • l] 3 



Tr 



(2^) 3 [ W 2 + ^l 2 ] 2 [( w + fi„) 2 + 4(1 + q) 2 

at g 2 r2 2 



16 (4g 2 + 2 )l 



(3.20) 



Next, we evaluate 



Da.a + D. 



3 + ^8.4 



2^4 iu — — + 2i( u> 



(3.21) 



In Eq. (3.21), the log-divergent constants A' and _B' fol- 
low from Eq. (3.5): 



A' = B' = -G s hxA. 



(3.22) 






FIG. 9: Corrections to the gauge (Coulomb) field self-energy 

too(i{j;,g;}). 



Diagrams D 8 . 5 and D 8 . 6 can be evaluated by simply 
multiplying the divergent vertex correction [Eqs. (3.12) 
and (3.13)] (evaluated at the left and right sides of the 
two- loop graphs, as shown in Fig. 10) by a factor of the 
non-divergent polarization bubble, given by the N = 00 
photon self-energy [Eq. (2.20)]. We obtain 



D 8 . 5 ~27 ? w 2 / 1 (M;)lnA 



2 



D 



8.6 



-2wg s In A 



v F q- 



(3.23a) 
(3.23b) 



Summing the results of Eqs. (3.18), (3.21), and (3.23) 
with the inverse N = 00 gauge propagator [Eq. (2.22)], 
and using Eqs. (3.19) and (3.22), we obtain the bare ver- 
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11.3 



D 



11.5 



D 



11.2 



D 



11.4 



D 



11.6 



Diagrams Dix.j-Dii.4 describe the autorcnormaliza- 
tion of the disorder at one loop. One finds 19 



Dai + D11.2 + D u .3 + D11.4 

2g u Q s + &g m gA + ^9v9A3 

. + 2g m (g~u + 9v). 

2irv 2 F (a a k 3 ® a a k 3 ) [2g 2 m + 8g 2 A + 2g u g v ] In A 



--2nv% (1 ® l) 



In A 
In A 



+ 2nv F (a 6 kp <g> a hp) 
+ 2irv 2 F (a 3 k 3 ® a 3 k 3 ) 



2gm(gv - g u ) 

+ 4g A (g u + g v ) 
2g v (2g m - g u - g 

+ 4gA3(gu + g v ) 
+ 8g A (g m - gv) 



In A 



In A, 
(3.27) 



FIG. 11: Renormalization of the disorder strength at 
O {j?,g 2 ), where g € {gZ, G^}. 



where we have re-introduced the following notation: 
barred Greek indices run over the "spatial" Pauli matrix 
components 



tex function 

*r<°) = - 



q - to 



v F q 



VFq 



+ 



VFq 2 n 2 n 



2 

■qw 2 f 3 (w) - wQ s 



In A 



-Vw 2 f 3 {w) + wQ 6 



In A. 

(3.24) 



Since there is no correction proportional to |q|, we see 
immediately that 



d In Z a 
dlnA 



0, 



(3.25) 



as assumed above in Eq. (3.15). The renormalization 
condition is then 



din A 



*t<°> = o. 



(3.26) 



Eq. (3.26) recovers the previously-derived results, 
Eqs. (3.10) and (3.17). 



D. Disorder strength renormalization 

Finally, we compute the renormalization of the disor- 
der strength, using the diagrams pictured in Fig. 11-13. 
These diagrams correct the bare disorder vertex Dx 4, pic- 
tured in Fig. 1, with the tree- level amplitude given in 
Table I. 



Me {1,2}, 



(3.28) 



to be distinguished from unbarred Greek indices which 
run over all three Pauli matrix components, e.g. 



a, /3e {1,2,3}. 



(3.29) 



Diagrams Dn.5 and Du.g give identical contributions, 
involving the dressing of the disorder vertex by the 
Coulomb interactions at O (1/N). Their evaluation pro- 
ceeds similarly to that of the electronic self-energy (de- 
tailed in Appendix A). We obtain 



D11.5 + Dn.6 

=2ng u v 2 F (1 ® i) [-27710/1(10)] In A 

+ 2irg A v 2 F (o & kp ® cr & kp) [2r)wf 2 (w)} In A 
+ 2itg A3 v 2 F (a & k 3 ® a & k 3 ) [2rjwf 2 (w)} In A 



+ 2-Kg m v 2 F (a 3 kp 
+ 2ng v v F (a 3 k 3 £ 



' a 3 ki 



a 3 k 3 ) 



2r)wf 2 (w) 
+ 2r]wf 3 (w) 
2r]wf 2 (w) 
+ 2riwf 3 (w) 



In A 

In A. 

(3.30) 



The functions fi(w)-f 3 (w) were defined by Eqs. (3.3) 
and (3.4). 

The remaining diagrams pictured in Figs. 12 and 13 
involve 3-electron loops. We consider first the ring dia- 
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and 






D12.7 



FIG. 12: Renormalization of the disorder strength at 
O {jf,g 2 ), where g € {gZ, Q^}. 



grams in Fig. 12, beginning with D12.1: 



dtod 2 ld 2 Y [a"k v v F &-Y] ab \ 



(2^) 5 v 2 F \' 2 |l'l 

' a^k y [iuj + v F a ■ (1-1')]' 



Tr 



x [iuj + v F a ■ l]a A K; 7 
t x [iuj + v F a ■ 1] 
[lo 2 + v 2 f \ 2 ] 2 [lo 2 + v 2 F {\-\'f] 



(3.31) 



Simplifying, we obtain 



D12.1 =(^) 2 g;g v x v 2 F (^k,a a ) atb (& x k,) c . d 



2 5 iw 



1 + w 



where 



K?' = 



J /fi (l-l')0W a 



, (3.32) 



dujd 2 \d 2 V 

(2TT) 5 |l'|3[ W 2+l 2 ]2[ LJ 2 + ( 1 _ 1 /)2]' 

12 lnA ^^, (3.33) 



2 10 2tt 



K 



dujd 2 \d 2 V 



r«(i-r) / W 



1 In A 



)10 



(2tt) 5 |1'|3[ U ;2 + 1 2 ]2[ w 2 + ( 1 _ 1 /)2] 



2- 



(3.34) 



Eqs. (3.33) and (3.34) may be obtained by a) perform- 
ing the ultraviolet-convergent 3-electron loop momentum 
integral over 1 via Fcynman parameters, b) evaluating 
the frequency integration over the real line, and c) com- 
pleting the final ultraviolet-divergent momentum integral 
over 1 . 

Combining Eqs. (3.32)-(3.34), and carefully summing 
barred (unbarred) indices over {1, 2} ({1, 2, 3}), we finally 
obtain 



D 



In A 



12.1 



' 2tt (a 3 kp ® a 3 kp) [-Ag m g A ] 
y + 2ir (& 3 k 3 ® a 3 k 3 ) [-4g v g A3 ] t 



(3.35) 

The remaining diagrams in Fig. 12 give identical con- 
tributions. Next, we consider the diagrams shown in 
Fig. 13. Each of the diagrams D 13 i-D 13 4 represents an 
autorenormalization of the bare scalar potential disorder 
parameter g u , mediated by the Coulomb interaction w. 
Due to the 3-ring cancelation property encapsulated by 
Fig. 7, the sum of these diagrams gives (exactly) zero. 
The bare vertex function is 



iTf =2tt 9u v 2 f (i ® 1) 



+ 2Trg»v F (a»k u )< 
+ D u + D12, 



(3.36) 



where D m , to g {11, 12}, denotes the sum of all diagrams 
in Fig. to. 

The renormalization condition is 



dlnZ, 



din A 



din A 



*rg> = 0, 



(3.37) 




'13.1 




'13.3 




'13.4 



FIG. 13: Renormalization of the disorder strength at 
(§,g 2 ), where 



13 



which leads to the rcnormalization group flow equations 



A. One-loop flow equations 



dg u 
dl 

dgA 
dl 

dgAz 
dl 

dg m 
dl 



dg v 
dl 



~-2g u Q 8 -riwf 3 {w) + Sg m g A + ig v g A3 , (3.38a) 
-8g A gA3 + 2g m (g u + g v ), (3.38b) 
--2g 2 m + 8g A + 2g u g v , (3.38c) 



-ig m {gv - 9u) + 4gA(g u + g v ) + 2g m riwf 3 (w 
16w 



1 + w 

) + 4g A3 (g u * gv) 



+ 8g A (g m - g v ) + 2g v rjwf 3 (w) 



16w 
1 + w 



(3.38d) 



gvgA3, 

(3.38e) 



where we have used Eqs. (3.4), Eqs. (3.9), and (3.17). [Q s 
was defined by Eq. (3.6).] 



IV. ANALYSIS OF THE RENORMALIZATION 
GROUP FLOWS 



In this section, we use the results of our one-loop renor- 
malization group (RG) study, detailed in Sec. Ill, to at- 
tempt to understand the topology of graphene's phase 
diagram in disorder-interaction coupling strength space. 
Our primary goal is to determine what kind of low-energy 
theory one should use to understand the macroscopic 
electronic properties of a graphene sheet, using the RG 
as our guide. If the RG reveals a theoretically tractable 
(e.g., pcrturbatively accessible) phase or critical point, 
then we can use that framework to make predictions that 
can be tested against experiment. As we will see, the sit- 
uation appears to be more complex (at least to lowest 
order in 1/N), and the ultimate low-energy physics (at 
zero temperature) likely requires a description very differ- 
ent from the weakly-perturbed, massless Dirac electron 
picture used as our starting point here. 

Nevertheless, the RG does allow us to identify several 
possible scaling regimes; our most interesting result con- 
cerns the apparent robustness of one such scaling regime, 
at least as viewed from the vantage point of the weakly- 
disordered Dirac electron theory. This (crossover) regime 
turns out to be dominated by the non-Abelian vector po- 
tential disorder, characterized by the parameters g A and 
g A3 [Eqs. (2.9) and (2.11)] which appear, e.g., in the de- 
scription of clastic lattice deformations ("ripples"), 5,11-13 
as well as topological defects in the graphene lattice 
(Sec. II A 2 and Refs. 20,36). This principal result is 
demonstrated via numerical integration of the RG flows, 
and is discussed in subsection IV E, below. Potential ex- 
perimental manifestations of our results are discussed in 
Sec. V. 



Using rcnormalized perturbation theory, we obtained 
in the previous section the one-loop flow equations for 
the coupling constants of the large- N graphene field the- 
ory defined by Eq. (2.19). Six coupling strengths de- 
fine that model: the Fermi velocity vf, the dimension- 
less Coulomb interaction strength w [Eqs. (2.23) and 
(2.24)], the screened scalar potential disorder strength g u 
[Eq. (2.26)], as well as g A , g A3 , g m , and g v [Eqs. (2.8)- 
(2.11)] . The latter four disorder parameters were encoded 
in the disorder metric Q", defined by Eq. (2.14). From 
Eqs. (3.10), (3.17), and (3.38) obtained in Sec. Ill, the 
flow equations to the lowest nontrivial order in the small 
parameters 1/N and {gu,Q^} are given by 



din vf 

dl 
dlnw 

dl 
dg u 
dl 



z-l-Q a + r)wf 3 (w), 

Qs - T}wfo(w), 

1 + w 



Qs - rjwf 3 (w) 



(4.1a) 

(4.1b) 
8g m gA + 4g v gA3 



{l + wy 



—jt = 8gAgA3 + 2g m (gu + g v ), 



dgA 
dl 

dgA3 
dl 

dg m 
dl 



dg v 
dl 



25 , m + 8g| + 23 u 3„, 



(4.1c) 
(4.1d) 

(4.1e) 



2gm(gv - g u ) + ^gA{g u + g v ) + 2g m vwf3(w) 

- 16 " gmgA , (4.1f) 
1 + w 

16w g v gA3 



+SgA(g m - g v ) + 2g v r]wf 3 (w) - 



In these equations, 



and 



1A + 2g A 3 + 2g„ 



gv, 



1 + w 
(4-lg) 



(4.2) 



ttN 



(4.3) 



In Eq. (4.1a), z = z(l) is the (as yet unspecified) dy- 
namic critical exponent. The function f 3 (w) appearing 
in Eq. (4.1) was defined by Eq. (3.4), repeated below for 
convenience: 



f 3 (w) 



-—{- 

2w 2 [it 



w 



arccos(ui) 
VI - w 2 



(4.4) 



valid for w < 1. The analytical continuation of f 3 (w) for 
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FIG. 14: Plot oiwf 3 (w). 
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w > 1 is given by 



2«r 



2w 

- 1 + — 

7T 

2 



= hi(w + \/ w 2 — lj 



(4.5) 



In the small and large w limits, f 3 (w) has the behaviors 

. . 1 2w 3w 2 n ^ _ 1 

f 3 ( w ) ~ -- — + —— + ... 0<w<l, 



— + 0[^), w>l. (4.6) 



7TW 



((!- 



The function wf 3 (w) is plotted in Fig. 14. 

Eq. (4.1c) follows from the definition of the screened 
scalar potential disorder strength g u [Eq. (2.26)], and 
the flow equation (3.38a) for the bare scalar potential 
strength. We emphasize that it is g u that appears in 
physical quantities in the large-iV theory, as demon- 
strated in Sec. V; the renormalization group must there- 
fore be formulated in terms of g u , the effective scalar 
potential strength. 



B. Termination of the RG flow for finite 
temperatures and chemical potential 

The derived RG equations (4.1) are limited by the 
regime of the weak disorder strength and usually become 
non-usable when I — ► oo. There are, however, physi- 
cal situations when those equations are adequate to de- 
termine the transport properties of the system. Those 
are regimes of either sufficiently large doping (charac- 
terized by the chemical potential //), or of sufficiently 
high temperature T. At such an energy scale, the log- 
arithmic renormalizations must be stopped, and kinetic 
coefficients computed, as in Sec. V. 

To use Eqs. (4.1) we have to relate the scale Z* at which 
the RG is stopped to the energy scale 



max (\fx\,T) . 



(4.7) 



To achieve this goal, we write a scaling equation 
dine* . . 

: -Z(l*), 



dL 



(4.8) 



and find the scale-dependent dynamic critical exponent 
by setting Eq. (4.1a) equal to zero. 34 ' 51 It gives us an 
implicit function Z*(e») 



=- f'dlz(l); 
v F Jo 

z(l) = l + g s (l)-r,w(l)f s [w(l)}, 



(4.9) 



where is the Fermi velocity defined on the spatial 
scale of the order of lattice constant a. 



C. Restricted flow equations: SU(2) valley space 
rotation symmetry restored 

A useful subset of the RG flow provided by Eq. (4.1) is 
to consider the system with SU(2) valley (k) space rota- 
tional symmetry restored on average. Such a restriction 
constrains 



9 A = 9 A3, 
9m 9v • 



(4.10) 



[See Eqs. (2.9) and (2.11).] Imposing Eq. (4.10), the RG 
equations reduce to 



din vf 

dl 
din w 

dl 
dg u 
dl 



dgA 
dl 

dg m 
dl 



z-1 

{(hi + &9A + 3.g. 
2g, 



1 + VJ 

Y2g m gA 



Oh, + ®9A + 3g m ) + r)wf 3 (w)(4.11&) 
-T)wf 3 (w), (4.11b) 

(4.11c) 
(4.11d) 



[gu + 6g/ 



(l + w) 2 ' 
8g 2 A + 2g m (g u + g m ), 

2g m {g.. 



9u) + 4g A 9u + 4 



1 - 3w 
l + w 



+2g m r]wf 3 (w). 



g m gA 
(4.11e) 



As discussed below in Sec. IV E, statistical SU(2) Fermi 
space (k) rotational symmetry is generically restored un- 
der integration of the full flow equations, Eq. (4.1); we 
will therefore focus upon the restricted flow equations 
(4.11) when we discuss our primary results in Sec. IV E. 

Note that Eq. (4. lid) implies that the only non-trivial 
fixed point structure (a fixed line) occurs in these re- 
stricted equations when g A = g m = (since the disor- 
der parameters, being variances [Eq. (2.11)], cannot take 
negative values). 
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D. Graphene with non-generic disorder 



2. Particle-hole symmetry and the chiral class BDI 



Before analyzing graphene subject to generic (short- 
range correlated), time- reversal invariant disorder, we 
discuss several simpler limiting cases that arise upon fine- 
tuning the disorder distribution defined by Eqs. (2.9) and 
(2.11). To do so, we specialize the flow Eq. (4.1) to the 
quantum disorder universality classes discussed in sub- 
sections IIA1-IIA3. Although it is typically difficult 
to fine-tune the disorder profile experimentally, some of 
the limiting cases discussed below may dominate vari- 
ous scaling or crossover regimes in graphene. (See also 
Sec. V.) 



1. Scalar potential disorder and the sympletic class All 

If remote charged impurities provide the principal 
scattering mechanism in the experiments discussed in 
Refs. 1-4, one would expect fluctuations in the poten- 
tial u(r) to dominate over the other disorder types. As 
a zeroth order approximation, we may neglect the other 
disorder potentials altogether, setting gA = Qaz = 9m = 
g v = [Eqs. (2.9) and (2.11)]. As discussed in Sec. II A 1, 
this is equivalent to enforcing SU(2) valley space rota- 
tional symmetry [invariance of Eq. (2.9) under K-space 
rotations] in every realization of the static disorder. This 
theory satisfies the effective TRI* condition, defined by 
Eq. (2.7), and belongs to the symplectic ("spin-orbit") 
ordinary metal class AIL 19-24 The flow equations are 



dim 



=g u -nwf 3 (w), 
1g u 



dg u 
dl ~l + w 



9u - r)wf 3 (w) 



(4.12a) 
(4.12b) 



These equations possess an unstable fixed line for g u 
r\wjz(w)\ moreover, 



g u (l)=Cw 2 (l), 



(4.13) 



where C — g u (0)/w 2 (0), so that the RG flows are 
parabolic trajectories in the g u — w plane. [g u = g u {^ + 
w) 2 is the unscreened scalar potential disorder strength; 
see Eq. (2.26).] 

A similar result was found previously for the case of 
graphene at weak Coulomb interaction coupling; 32,33 in 
the case of the large- A generalization, a technical error 
in Rcf. 32 (see Sec. Ill A and Appendix A) led to the 
prediction that the line parameterized by Eq. (4.13) be- 
comes attractive at larger Coulomb interaction strengths. 
Instead, Eq. (4.12), coupled with the monotonic behavior 
of wf 3 (w), [Eqs. (4.4)-(4.6) and Fig. 14], shows that this 
line is repulsive for all w > 0. In other words, for the case 
of scalar potential disorder only, the physics in both the 
weak coupling and large-A pictures is qualitatively the 
same. The stability of the g u -w fixed line in the full 6- 
dimcnsional disorder-interaction coupling constant space 
is analyzed in Appendix B. 



As discussed below Eq. (2.11) in Sec. II A, different 
quantum disorder (random matrix) universality classes 
may be theoretically realized by enforcing invariance of 
the graphene system, in every realization of disorder, 
under different combinations of the transformations de- 
fined by Eqs. (2.4)-(2.7). If we enforce both particle- 
hole symmetry (PH) and time-reversal invariance (TRI) 
[Eqs. (2.4) and (2.5)] in every static disorder realization, 
consistent with, e.g., the presence of carbon atom vacan- 
cies (treated as scattering centers in the unitary, hard- 
scattering limit) and the absence of further-neighbor 
hopping, 52 then the system falls into the "chiral or- 
thogonal" class BDI. 20 ' 27-29,40 ' 41 The allowed disorder 
strengths are g m and 17^3, and the one-loop flow equa- 
tions are 



d In w 
dl 

dgA3 
dl 



~-2g A3 + 2g m -nwf 3 (w), 



--29L 



=2g m rjwf 3 (w). 



(4.14a) 
(4.14b) 
(4.14c) 



These equations possess no fixed point for g m ^ 0, and 
generically flow to both strong disorder (gA3,9m ^ 00 ) 
and strong interaction (w — > 00) coupling. In passing, 
we note that, in the absence of Coulomb interactions 
(w = 0), Eqs. (4.14) can be extended to all orders in 
g m and gAi, using conformal field theory methods. 28 One 
then finds that disordered graphene with strict particle- 
hole symmetry [Eq. (2.4)] and no long- or short-ranged 
interparticle interactions 34 possesses a critical, delocal- 
ized phase, indicative of the existence of extended single- 
particle states at zero energy, for arbitrary strength dis- 
order. 

The theory with g m = does possess an attractive 
fixed line in the gA 3 -w plane. This line has been dis- 
cussed before; 32,33 it is unlikely to play an important role 
in real graphene physics, because it is highly unstable in 
other disorder directions, as shown in Appendix B. 



3. Vector potential disorder, ripples and topological defects: 
class CI 

If we enforce both PH* and TRI [Eqs. (2.6) and 
(2.5)], then the system falls into the CI quantum disor- 
der class. 20,40,41 Recall from Sec. II A that PH* describes 
an effective particle-hole transformation, different from 
the physical PH transformation defined by Eq. (2.4), 
the latter of which is inherited directly from the lattice 
model [Eq. (2.1)]. The allowed disorder strengths are gA 3 
and gA, i.e. the pure Dirac theory perturbed only by a 
quenched, SU(2) non-Abelian vector potential. 

As discussed in Sec. II A 2, the intravalley (Abelian) 
vector potential {A?} appears in the description of 
long wavelength "ripples" in the low-energy Dirac 



16 



theory. 5,11,12 Both the intravalley and intervalley {A^} 
vector potential components play a primary role in the 
field theoretic description of honeycomb lattice disloca- 
tions and disclinations. 36 
The RG flow equations are 



dim 



=ig A + 2g A3 -rjwf 3 (w), 



dgA „ 

—jj- =ogA9A3, 

dg 



.43 



dl 



A- 



(4.15a) 
(4.15b) 
(4.15c) 



These equations possess no fixed point for gA ^ 0, 
and generically flow to both strong disorder and inter- 
action coupling. In the absence of interactions, w = 0, 
the Dirac theory with only non-Abclian vector poten- 
tial disorder flows to strong coupling; the flow termi- 
nates at an exactly-solvable conformal field theory fixed 
point 18,26 [equivalent to the Sp(2n) principal chiral non- 
linear sigma model, augmented with a WZW term]. 40 
Like the non-interacting BDI class discussed above, the 
fixed point of the non-interacting CI class model (with 
both gA and gA3 nonzero) describes a critical, delocal- 
ized phase possessing extended wavefunctions near zero 
energy. 18,26 This disorder class (in the absence of in- 
teractions) was discussed specifically in the context of 
graphene in Ref. 20; these authors argued that the con- 
ductance of the non-interacting CI class model should be 
of order the conductance quantum, and independent of 
the disorder strength. 



E. Graphene: generic disorder 

At last, we turn to an analysis of our large- N graphene 
field theory, characterized by the full one-loop RG flow 
equations (4.1). In the absence of any fine-tuning of the 
disorder potentials, our graphene model possesses time- 
reversal invariance [Eq. (2.5)] and (physical) spin SU(2) 
rotational symmetry in every disorder realization. This 
system is in the standard orthogonal metal class AI. In 
the absence of interactions, all single particle wavefunc- 
tions are expected to be exponentially localized in 2D. 44 

The full flow Eqs. (4.1b)-(4.1g) possess several fixed 
line structures involving a single nonzero disorder pa- 
rameter and the Coulomb interaction strength w. The 
case of scalar potential disorder g u was discussed above, 
in Sec. IV D 1. Weak coupling analogs of this and other 
fixed lines were discussed in Refs. 32,33, while large-Af 
versions were obtained and analyzed in Ref. 32 [but see 
the discussion following Eq. (4.13), above]. We will not 
dwell upon these structures here, however, because all of 
them are exceedingly unstable in the full 6-dimensional 
disorder-interaction coupling constant space. (This fact 
is demonstrated explicitly in Appendix B.) 

The RG flow equations [Eq. (4.1)] possess no per- 
turbativcly accessible, isolated fixed points, other than 




FIG. 15: Set of 10 RG flow trajectories obtained via numerical 
integration for different initial Coulomb interaction strengths 
w. Here, we study the case with SU(2) valley space symme- 
try restored on average, using Eq. (4.11). The initial condi- 
tions for the disorder coupling strengths are g~Z,(0) — 5a(0) = 
<?m (0) = 0.01 for all trajectories shown in the figure. We 
have set the expansion parameter n — 0.3 [see Eq. (4.3)]; the 
shape of the trajectories is only weakly dependent upon n 
(within our one-loop approximation). Lower-case Latin let- 
ters label trajectories with different initial Coulomb interac- 
tion strengths, given by w = (a) 0.010, (b) 0.079, (c) 0.13, (d) 
0.16, (e) 0.20, (f) 0.25, (g) 0.32, (h) 0.40, (i) 0.63 (j) 1.0. We 
do not exhibit the RG evolution of the random mass parame- 
ter g m [see Eqs. (2.9), (2.11), (4.11)], because we find that its 
behavior is always subleading compared to the scalar (g~Z) or 
vector (g A ) potential parameters. 



the (trivial) non-interacting, clean (not disordered) Dirac 
fixed point, defined by the condition 



9u = 9 A = 9 A3 = 9m = 9 



= 9 m = 9v = 0. 



(4.16) 



In order to gleam information about graphene's low- 
energy physics (to order 1/N in our large- N calculational 
scheme), it is necessary to integrate the flow equations, 
and observe the evolution of the flow trajectories as a 
function of initial conditions in coupling constant space. 
In doing so, we seek to address two crucial questions that 
so far remain unanswered: 

1. Over what range of initial conditions, if any, does 
the theory flow back to the clean, non-interacting 
Dirac fixed point [Eq. (4.16)]. In other words, can 
we find a critical manifold of codimension 0, with 
the Dirac fixed point as its sink? 

2. If the clean, non-interacting Dirac fixed point is 
found to be unstable, where does the system flow 
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FIG. 16: Same RG trajectories shown in Fig. 15, but pro- 
jected into the g~Z — gA disorder plane. 

to? What is the correct low-energy theory that 
we should use to understand disordered, interacting 
graphene? 

We can easily answer the first question, but we will be 
forced to speculate regarding the second. 

We have numerically integrated Eqs. (4.1b)-(4.1g) over 
a wide range of disorder and interaction coupling strength 
initial conditions. In our one-loop approximation, the 
resulting flow topology is only weakly dependent upon 
the expansion parameter rj, defined by Eq. (4.3). Wc 
find that the non-disordered, non-interacting Dirac fixed 
point, located by Eq. (4.16), as well as the entirety of 
the adjoining clean, but Coulomb-interacting line, < 
w < oo and g u = 9A = 9 A3 = 9m = 9v — 0, is unstable 
in the presence of arbitrarily weak, but generic disor- 
der [i.e. all five disorder strength parameters in Eq. (4.1) 
nonzero, but arbitrarily small]. We stress that the en- 
tirety of the clean, interacting line w > is perturba- 
tively accessible in the N — > oo limit, and flows back to 
the non-interacting Dirac fixed point in the absence of 
disorder. 31,35 

In the advent of nonzero quenched randomness, the 
flow to strong coupling exhibits two key features: a) 
SU(2) valley space rotational symmetry is generically 
restored on average upon integration of Eq. (4.1) for 
general initial conditions, i.e. the differences (\gA — 
9A-i\)/{9A + 9as) and (\g m - g v \)/(g m + 9v) asymptote 
to zero as the solutions to the RG equations reach the 
limit of their validity. From here on in, we will therefore 
restrict our discussion to the flows of the valley space 
SU(2)-symmetric system, described by Eq. (4.11). b) 
The direction of the diverging RG flows depends strongly 
upon the initial strength of the Coulomb interaction pa- 
rameter w. For weak (or vanishing) Coulomb interactions 
W <C 1, the flow is dominated by the divergence of the 



screened scalar potential disorder parameter g u [defined 
by Eq. (2.26)]. The ratio of the other disorder strengths 
9a and g m to g u asymptotes toward zero as the flow 
begins to diverge. For stronger values of the Coulomb 
interaction strength w > 0.1 (for the choice rj = 0.3), 
we observe a "crossover" in the flow direction. In the 
limit of large Coulomb strength w > 1, we find that the 
flow becomes dominated by the divergence of the SU(2) 
non-Abelian vector potential disorder parameter gA\ in 
this regime, the ratio of the other disorder strengths g u 
and g m to gA asymptotes to zero as the flows leave the 
perturbatively accessible regime. 

We demonstrate the picture described above with a 
selection of RG flow trajectory plots. Figs. 15-18 de- 
pict two sets of flow trajectories obtained via numeri- 
cal integration of Eqs. (4.11b)-(4.11e) for different initial 
Coulomb interaction strengths w. Figs. 15 and 17 depict 
projections of these flows in the 3D coupling constant 
subspace {g u , gA, w). Figs. 16 and 18 respectively exhibit 
2D projections of the same flows represented in Figs. 15 
and 17 in the disorder plane (g u ,gA)- We choose not to 
exhibit the RG evolution of the random mass parame- 
ter g m [Eqs. (2.9), (2.11), and (4.11c)] in these figures, 
because we find that g m always plays a subleading role 
relative to either the scalar g u or non-Abelian vector gA 
potential disorder strengths. 

In each of Figs. 15-18, all trajectories share a given set 
of initial disorder strengths (described in the figure cap- 

log 10 w 




FIG. 17: As in Fig. 15 but with the initial conditions for the 
disorder coupling strengths g~Z(0) = <?m(0) = 0.1, gA = 0.01. 
Lower-case Latin letters label trajectories with different initial 
Coulomb interaction strengths, given by w — (a) 0.010, (b) 
0.25, (c) 0.40, (d) 0.50, (e) 0.63, (f) 0.79, (g) 1.0, (h) 1.3, (i) 
1.6, (j) 2.5, (k) 6.3. 
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FIG. 18: Same RG trajectories shown in Fig. 17, but pro- 
jected into the gZ — gA disorder plane. 



tions). Different trajectories, distinguished by lowercase 
Latin letters, correspond to different initial Coulomb in- 
teraction strengths w(0). Trajectories with successive la- 
bels (a), (b), etc. carry successively larger initial values of 
w; numerical values are stated in the captions of Figs. 15 
and 17. Although very weak initial Coulomb interaction 
strengths w < 0.1 lead to trajectories that are dominated 
by large scalar potential disorder fluctuations g u — > oo, 
stronger initial interaction strengths bend the trajecto- 
ries away from the g u -w plane, toward the particle-hole 
symmetric gA~w plane. [The effective particle-hole sym- 
metry is defined by the PH* transformation given by 
Eq. (2.6)]. 

These figures encapsulate our primary result: non- 
Abelian vector potential disorder dominates the flow to- 
ward strong coupling (gA, w — > oo) for moderate to 
strong Coulomb interaction strengths. As stated above, 
this result is essentially independent of our expansion pa- 
rameter ?7, defined by Eq. (4.3), for our one-loop flow 
equations (4.1) and (4.11), obtained to lowest nontriv- 
ial order in 1/N and in the disorder strengths g u , gA, 
g m . Physically, strong scalar potential disorder fluctu- 
ations would favor the accumulation of electrons and 
holes in spatially segregated "puddles," locally violat- 
ing charge neutrality. The advent of a second type of 
disorder that manifestly preserves a type of particle-hole 
symmetry [Eq. (2.6)] allows the system to flow toward a 
strongly disordered, interacting regime, while everywhere 
preserving electroneutrality. Equivalently, Coulomb in- 
terparticle interactions lead to screening effects, which 
parametrically curtail the scalar potential disorder fluc- 
tuations; the non-Abelian vector potential is not affected 
by the interactions because it represents charge-neutral 
randomness (such as the vector potential component of 
"ripples" ). 
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FIG. 19: Schematic phase diagram for disordered, interact- 
ing graphene, projected into the disorder (gZ, <7a) interaction 
(w) coupling constant space, based in part upon the pertur- 
bative RG calculations performed at weak coupling, 19 ' 33 and 
within the large- N approximation (as detailed in this paper) . 
The labeled small dots represent "known" phases, while the 
large black globes represent "non-linear sigma model (NLcrM) 
arenas," wherein the system is described by a well-known (if 
poorly understood) model, which may possess multiple phases 
in d = 2 spatial dimensions. 44,51 Point (a) represents non- 
interacting, clean (not disordered) graphene, as located by 
Eq. (4.16). This theory is a sink for the clean, Coulomb 
interacting line (the vertical axis), as well as a portion of 
the g~Z — w plane. It is unstable to generic disorder per- 
turbations. The dashed line in the g~Z — w is the unstable 
fixed line located by Eqs. (4.12a) and (4.12b). The point 
(b) is the non-interacting, non-Abelian SU(2) vector poten- 
tial disorder CFT. 18 ' 26 Both points (a) and (b) are embedded 
in the particle-hole symmetric plane g~Z = g m = 0, repre- 
sented by the shaded plane in the figure. Point (c) repre- 
sents the non-interacting Anderson insulator for the orthogo- 
nal metal class (AI) in 2D. Globe (d) should be described by 
the symplectic/spin-orbit class (All), non-interacting NLcrM, 
possibly augmented by a topological term. 19-24 In princi- 
ple, both diffusive metallic and Anderson insulating phases 
are possible. 44 Globe (e) represents the symplectic class 
Finkel'stein non-linear sigma model (FNLcrM) for interacting 
electrons; 51,53,54 globe (f) represents the orthogonal class (AI) 
FNLo-M. 51,55,56 The numerical RG flows exhibited in Figs. 15 
and 17 should be "superimposed" upon this figure. 

To close this section, we direct the reader's attention 
to Fig. 19. In this figure, we exhibit a putative phase 
diagram for graphene, projected into the 3D (g~ u ,gA,w) 
disorder-interaction coupling strength space, based upon 
our work here, as well as previous results 18,20,51,55,56 re- 
garding the zero temperature ground state physics of 2D 
disordered and/or interacting electronic systems. The 
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most important features of this figure are the points la- 
beled (a)-(c) and the "globes" labeled (d)-(f). Points 
(a)-(c) label known phases of the graphene model in the 
absence of interactions: (a) is the non-disordered, non- 
interacting Dirac theory, (b) denotes the critical, delo- 
calized phase that occurs at strong non-Abelian vector 
potential ("chiral") disorder, described analytically by 
the class CI conformal field theory 18,26 As discussed in 
Sees. II A 2 and IV D 3, the class CI fixed point is ex- 
pected to possess a non-zero conductance at zero tem- 
perature that is independent of the disorder strength. 20 
Finally, (c) denotes the Anderson insulating phase for 
the orthogonal normal metal quantum disorder class. 
Globe (d) represents a "theory arena" corresponding to 
the symplectic (spin-orbit) normal metal class for non- 
interacting electrons, described by the symplectic non- 
linear sigma model (possibly modified by a topological 
term). 21-24 The symplectic class may possess both metal- 
lic and insulating states in 2D. 44 

Globes (e) and (f) respectively describe the sym- 
plectic and orthogonal normal metal quantum disorder 
classes, augmented with long-range Coulomb interpar- 
ticle interactions. 51,53-56 These globes describe physics 
on temperature scales less than the elastic scattering 
rate, i.e. T < 1/t e ;, with T e i the elastic transport life- 
time due to impurity scattering. The low-energy field 
theory description takes the form of a modified non- 
linear sigma model that incorporates electron-electron 
interactions. 51,55 The Drude conductivity computed in 
the regime T ~ T e i, discussed in the next section, enters 
through the bare diffusion constant of the sigma model. 

Arrows in Fig. 19 schematically denote RG flows whose 
topology has either been derived in this paper, previous 
work, 51 or may be inferred by general principles. The 
dashed line residing in the g u -w plane represents the re- 
pulsive fixed line described in Sec. IV D 1; the section of 
the g u -w plane between this line and the Coulomb in- 
teraction (to) axis drains into the clean, non-interacting 
Dirac fixed point (a). This section is, however, highly 
unstable in the other disorder (i.e. qa) directions. The 
reader should imagine "superimposing" the flow trajec- 
tories depicted in Figs. 15 and 17 upon this figure. 



V. PHYSICAL RESULTS AND DISCUSSION 

In this section, we articulate a number of predictions 
for the scaling behavior of the dc electrical conductiv- 
ity, as well as the thermal transport, for graphene in its 
Drudc/Boltzmann transport regime. The predictions are 
extracted from the large-TV, weak disorder renormaliza- 
tion group (RG) analysis performed above. We presup- 
pose here that the temperature is sufficiently high so that 
quantum interference corrections to the conductivity due 
to electronic diffusion (weak localization and Altshuler- 
Aronov-type interaction effects) may be ignored. 44,57,58 

We restore Planck's constant via h = 2ir throughout 
this section. 



Disorder vs. interaction limited transport 
coefficients 



The kinetic transport coefficients in graphene are ex- 
plicit functions of the effective disorder and interaction 
strengths, depending upon both the elastic transport life- 
time r e i due to impurity scattering, and upon the inelas- 
tic transport lifetime Tin due to electron-electron colli- 
sions. Estimates for these are 



— oc Q T ma,x(\ fj,\,T), 

Tel 



(5.1a) 



h l t 2 

oc — min(w 2 , 1) min(T, — ), (5.1b) 

Tin A 

where /i is the chemical potential, T is the temperature, 
and Qj is a certain combination of the disorder parame- 
ters in Eq. (2.11)-see Eqs. (5.3) and (5.4), below. [Recall 
that to tx r s , Eq. (2.24).] 

We consider the two limits with large and small T e i /n n 
in turn: 



1. r e i > Ti n : interaction- limited 

Real, inelastic electron-electron collisions play an im- 
portant role in determining the kinetic coefficients for 
T e i > Ti n . In particular, at exactly zero doping \i = 0, 
the composite electron-hole fluid is electrically neutral. 
Here, electron- hole collisions are enough to set a finite dc 
conductivity Ode, which takes the form 



V 2 e 



2 C 2 



G dc OC 



h min(io 2 ,l) 



(5.2) 



For Coulomb interaction strengths to > 1, Eq. (5.2) pre- 
dicts a "minimum metallic conductivity" that is indepen- 
dent of both the interaction and disorder strengths (to a 
first approximation). The doping dependence of the con- 
ductivity is limited by the impurity scattering. 

Violations of Mott's formula [Eq. (5.6)] for the ther- 
mopower and of the Wicdemann-Franz law for the ther- 
mal conductivity are expected in the interaction-limited 
regime. Quantitative results based on the relativistic hy- 
drodynamics description of the electron-hole plasma in 
graphene will appear elsewhere. 59 



2. T e i < Tin'- disorder- limited 

When Tin due to electron-electron collisions exceeds 
T e i due to the impurities, the situation is simpler. To the 
first approximation, we can neglect contributions to the 
kinetic coefficients due to inelastic processes. The Drude 
conductivity in the ladder approximation then takes the 
form 



Ne 2 



Ode 



-E F Tel, 



(5.3) 
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where 59 

E F T e l 



1 



7T<5t 

Qt = flu + 8.9a + 4gA3 + Gg n 



3g v - 



(5.4a) 
(5.4b) 



Note that it is the screened scalar potential disorder 
parameter g u [Eq. (2.26)] that appears in Eq. (5.4). 
The temperature or chemical potential-dependence of 
the conductivity follows from combining the scaling pre- 
dictions of Sec. IV with Eq. (5.3); Mott's formula and 
Wiedemann- Franz apply for thermotransport. In the re- 
mainder of this section, we will use this strategy to iden- 
tify three scaling regimes for the transport coefficients in 
the disorder-limited case. 
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B. Scaling predictions for graphene transport 

The essential message is as follows: if the disorder 
distribution is sufficiently weak and short-range corre- 
lated, then the RG generically predicts corrections to 
the conductivity that are logarithmic in temperature. 
In particular, we have identified three scaling regimes 
which may be observed in (future) graphene experiments: 
(1) a "QED" regime, wherein the conductivity increases 
with decreasing temperature, as well as (2) intermediate 
and (3) "QCD" regimes, characterized by a conductivity 
that decreases monotonically with decreasing tempera- 
ture. Here, the descriptors "QED" and "QCD" refer 
to formal analogies between the disordered, interacting 
graphene theory in particular scaling regimes, and the 
theories of high energy quantum electro- and chromo- 
dynamics, respectively; we do not imply a mapping be- 
tween experimental phenomena in high energy particle 
and graphene physics. 

In this subsection, we set N = 4, appropriate to real 
graphene. To spare the reader from technical details, the 
derivation of the formulae presented below has been rele- 
gated to Appendix C. A sketch of our results is pictured 
in Fig. 20. 

1. "QED" regime: Scalar potential disorder and flow 
toward weak coupling 

The bare Coulomb interaction strength r s = e 2 /ehvp 
is of order unity in the previous substrate-supported 
experiments, 1-4 and should exceed two for suspended 
films. 11 As discussed in Sees. I and IV, an r s ~ 1 rep- 
resents a relatively strong coupling regime for the mass- 
less Dirac electrons in undoped graphene. If the disor- 
der effects are weak relative to the Coulomb interaction 
strength, then the RG predicts the existence of a "QED" 
scaling regime, wherein disorder and electron correlation 
effects grow ever weaker on larger length or lower en- 
ergy scales. This scaling toward weak coupling manifests 
itself as a logarithmic-squared increase in the dc conduc- 
tivity with decreasing temperature [Eq. (5.5), below]. A 



FIG. 20: Sketch of possible logarithmic scaling regimes of 
the dc conductivity <7dc vs. temperature in undoped graphene. 
The labels (l)-(3) identify the qualitative scaling behaviors 
pictured here with the three regimes respectively discussed in 
subsections V B 1-V B 3 in the text. Non-monotonic behavior, 
i.e. the peak separating the curve segments labeled (1) and 
(3) in the figure, can occur if a sample microscopically domi- 
nated by strong Coulomb interactions and weak scalar poten- 
tial disorder ("QED" regime), distinguished by a conductivity 
that increases with decreasing temperature [Eq. (5.5)], crosses 
over at lower temperatures to the strong coupling ("QCD") 
regime dominated by non-Abelian vector potential disorder, 
characterized by a decreasing conductivity [Eq. (5.10)]. An 
intermediate scaling regime (2) is also possible, characterized 
by a monotonically decreasing conductivity [Eq. (5.8)]. These 
predictions require that r ei < Ti„; see text. 

further requirement to reach the "QED" scaling regime 
is that scalar potential disorder must provide the domi- 
nant scattering mechanism for the Dirac electrons. (Con- 
sult Sees. IIA1-IIA3, IVD 1, and the last paragraph of 
Sec. IV E.) 

In this "QED" regime, we find that (see Appendix C 
for details) 

where Tr 3> T is an arbitrary reference temperature, and 
9~u <C r s ~ 1 is the microscopic strength of the (screened) 
scalar potential disorder [Eqs. (2.9)-(2.11), (2.26)]. 

The behavior described by Eq. (5.5) cannot per- 
sist indefinitely, because the presence of other disor- 
der parameters will eventually induce a crossover to the 
"QCD" regime described below in subsection VB 3. This 
crossover would be distinguished by a non-monotonicity 
of the conductivity, as shown in Fig. 20. We do not pro- 
vide here an estimate for the non-universal crossover tem- 
perature T c , which is difficult to extract quantitatively 
because it depends upon the details of the full micro- 
scopic disorder distribution. 60 

We also discuss thermotransport. We neglect interac- 
tion contributions or corrections to the energy current of 
the Dirac electrons; as a result, we find that the thermal 
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conductivity is slaved to Eq. (5.5) via the Wiedemann- 
Franz law. Within this same approximation, the ther- 
mopower S is given by Mott's formula 61 (fcs = 1): 
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(5.6) 



In this expression, (Jdc(e) is the scaling form of the con- 
ductivity, with temperature replaced by e (i.e. the chem- 
ical potential). Using Eq. (5.5), we obtain 
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for \/j,\ > T, with fj,^ > /i. 



(5.7) 



2. Intermediate regime: Scalar potential disorder and flow 
toward strong coupling 

We next consider the case where the strength of the 
Coulomb interaction r s is made weaker than that of the 
screened scalar potential disorder fluctuations [g u > w ~ 
r s }. We also require that scalar potential fluctuations 
provide the dominant source of scattering for the Dirac 
electrons, while not becoming so strong so as to invalidate 
our weak disorder assumption (g u < 1). Such a scenario 
might prove difficult to realize experimentally, but could 
be achieved in principle by screening the long-ranged 
Coulomb interactions with an external gate (brought into 
very close proximity to the graphene layer) , or by coating 
graphenc with a high-k dielectric material. 

The scaling behavior of the conductivity in this regime 
is given by 



<Tdc(T) 



4e2 
hn 



9u \T 



(5.8) 



for T R >T (Appendix C). 

The thcrmopower at non-zero chemical potential > 
T is given by 



3 e fi 



(5.9) 



Coulomb interaction parameter r s > 0.1. [See Sec. IV E 
for details.] In this regime, the graphene theory becomes 
formally analogous to 3 + 1-D quantum chromodynamics 
(QCD), in that the system becomes dominated by strong 
non-Abelian vector potential disorder and electron corre- 
lation effects (due to the Coulomb interactions). The es- 
sential difference between graphene and real QCD is that 
in the former (latter) case, it is a static, time-independent 
(dynamic, gluon-mcdiated) non-Abelian gauge potential 
that drives the system toward strong coupling at low en- 
ergies. 

Since r s > 1 in real graphene, we generically expect 
to reach the "QCD" regime at some finite temperature 
scale. We obtain the following formula for the asymptotic 
scaling behavior of the conductivity (Appendix C): 



a dc (T) 
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Tr 
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, (5.10) 



where gA is the (renormalized) strength of the non- 
Abelian vector potential disorder, and Tr > T. 
Eq. (5.10) is not expected to hold down to arbitrarily low 
temperatures. As we examine lower and lower energy 
scales, the effective vector potential disorder strength 
<?a(T) grows ever larger (Sec. IV E); as a result, the 
scaling behavior of Eq. (5.10) will mutate as gA evolves. 
[A more general scaling formula in the "QCD" regime 
is given by Eq. (C3) in Appendix C, which reduces to 
Eq. (5.10) in the limit g A < 1.] 

The thermopower at finite chemical potential \pb\ > T 
corresponding to Eq. (5.10) is given by 



S 



n 2 T'< 



PL 7T 2 - 2 



(5.11) 



Finally, we reiterate the following point: we have 
assumed that the inelastic transport lifetime due to 
electron-electron collisions, T in , exceeds the elastic life- 
time T e i due to the disorder. In graphene, the kinetic co- 
efficients are sensitive to the value of Ti n /r e /; we will dis- 
cuss the implications of the above three scaling regimes 
in the limit T;„/T e ; > 1 in a subsequent publication. 59 
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3. "QCD" regime: Non-Abelian vector potential disorder 
and asymptotic flow toward strong coupling 

Non-abelian vector potential disorder provides an- 
other potentially important scattering mechanism in 
graphene. 5,11 ' 12 This disorder type arises in the context 
of elastic deformations or "ripples" and in the context 
of topological lattice defects. 36 (See also Sees. II A 2 and 
IVD3.) 

Surprisingly, we have found in this paper that 
quenched vector potential disorder generically emerges 
as the dominant scattering mechanism at long length 
scales and low energies for Dirac electrons with a bare 
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APPENDIX A: COMPUTATION OF THE 
ELECTRONIC SELF-ENERGY TO O (■£). 

In this appendix, we derive the amplitude given by 
Eq. (3.1) for the diagram D5.1, pictured in Fig. 5. D5.1 is 
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the contribution to the electronic self-energy due to the 
Coulomb interactions at order 1/A. Using the Feynman 
rules from Sec. II B, one has 



D5.1 = i 



uP\ f dio d 2 \ i(u) n + u)) + vf& • (k + 1) 



NJ J (2tt)3 K+ W ) 2 + V 2,(k + 1) 2 
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|1| y/vj.P + LJ* + WV F \l\' 

(Al) 



Expanding in terms of the external frequency and mo- 
mentum, we obtain 
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where 



, dujd 2 l 1 y/l 2 + uj 2 \ 2 -uj 2 tk . 
h = I , n ^ — — , — ~ (A3a) 



dud 2 \ 1 VI 2 +w 2 ^ 2 

2\2 ' 



(2tt)3 Illv^ + ^ + tSllI («a + l 2 ) 
and where it is useful to define 



(A3b) 



I s = h + h 



dud 2 \ 1 VI 2 + oJ 2 



(2tt)3 |l| ^/W+Z? + w\\\ ( w 2 + 1 2y 



(A3c) 



After rescaling w — > , the momentum integrals are 
trivially done, so that 
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where we neglect ultaviolet-finite terms associated with 
the infrared region of the momentum integration. We 
then find that 
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where the functions /1-/3 were defined by Eqs. (3.3) and 
(3.4). Using Eq. (A5) in Eq. (A2), we recover Eq. (3.1). 



APPENDIX B: FIXED LINES AND LINEARIZED 
FLOW EQUATIONS 

The one- loop, large- A flow equations given by 
Eq. (4.1) in Sec. IV possess a number of critical fixed line 



structures, continuously connected to the non-disordered, 
non-interacting Dirac fixed point [Eq. (4.16)]. In this ap- 
pendix, we enumerate these structures, and determine 
their stability in the six-dimensional coupling strength 
space [c.f. Eq. (4.1)]. Each fixed line resides in a plane 
formed between the Coulomb interaction axis w and a 
single, nonzero disorder parameter. 



a. Scalar potential disorder 

As discussed in Sec. IV D 1, in the case of pure scalar 
potential disorder, i.e. only g u > among the five dis- 
order parameters appearing in Eq. (4.1), the large- A 
graphene model exhibits a repulsive fixed line in the g u - 
w plane. From Eq. (4.12), the fixed line is parameterized 
by the conditions 



9 A = gA3 = 9 m = 9v = 0, 



(Bla) 
(Bib) 



Since g u and 77 = 8/ 7r A arc expansion parameters for our 
large- A theory, the entirety of the fixed line described by 
Eq. (Bl) is perturbatively accessible, because the func- 
tion wjj,(w) asymptotes to a finite constant as w — > 00; 
see Eq. (4.6) and Fig. 14. 

Linearizing Eq. (4.1) about the line in Eq. (Bl), we ob- 
tain the following RG eigenvalues (as a function of w^): 
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The first two eigenvalues listed in Eq. (B2) respectively 
characterize RG flow trajectories along and through the 
fixed line within the g u — w plane; the latter two eigen- 
values listed in this equation characterize the RG flow in 
the transverse directions, and are each doubly degener- 
ate. The second and fourth eigenvalues listed in Eq. (B2) 
are positive for all finite > 0; both the g~ u (w^) 
fixed line, as well as the portion of the g u -w plane span- 
ning the region between this fixed line and the ui-axis are 
therefore highly unstable. 
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b. Abelian vector potential disorder 

Next, we consider the following fixed line in the gA3~w 
plane: 



9u = 9 A = 9n 



9v = 0, 
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Linearizing Eq. (4.1) along this line, we obtain the eigen- 
values 

{0, -p (0 ' f 3 (w (0) ), 2t)iii (0) /alw' '),^' ) / 3 (tc(°)), 
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The first two eigenvalues listed in Eq. (B5) respectively 
characterize RG flow trajectories along and through the 
fixed line within the gA3~w plane, where it serves as a 
(locally) attractive, critical extension of the clean Dirac 
fixed point [located by Eq. (4.16)]. The latter four eigen- 
values appearing in Eq. (B5) characterize the RG flow in 
the transverse directions, and show that this fixed line is 
highly unstable. 



c. Staggered potential disorder 

Finally, we consider a fixed line in the g v -w plane, 
parameterized by the conditions 



9u = 9 A = 9 A3 = 9rn = 0, 

ff (°)=^(°)/ 3 (^W). 
Linearizing Eq. (4.1) yields the RG eigenvalues 
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The first two eigenvalues listed in this equation re- 
spectively characterize RG flow trajectories along and 
through the fixed line within the g v -w plane; the latter 
four characterize the RG flow in the transverse directions. 
We observe that the g^\iv^) critical fixed line is also 
unstable for all finite w^°> > 0. 



APPENDIX C: DERIVATION OF PHYSICAL 
RESULTS 

In this appendix, we derive the scaling predictions for 
the dc conductivity stated in Eqs. (5.5), (5.8), and (5.10) 
of Sec. V. These predictions arise in the "QED," in- 
termediate, and "QCD" scaling regimes respectively, as 
defined in that section. 




FIG. 21: Schematic RG flow diagram, projected into the 3D 
coupling constant space formed from the screened scalar po- 
tential disorder strength (g~Z), the non- Abelian vector poten- 
tial disorder strength ((/a), and Coulomb interaction strength 
(w = 7rr s /2) [Eq. (2.24)]. This diagram should be com- 
pared to the flow plots, obtained via numerical integration 
of Eq. (4.11), depicted in Figs. 15 and 17 of Sec. IVE. In 
this figure, we have sketched characteristic flows, designated 
by the labels (l)-(3), corresponding to the "QED," interme- 
diate, and "QCD" scaling regimes discussed in Sec. V. The 
dashed curve represents the unstable fixed line discussed be- 
low Eq. (4.12) in Sec. IV D 1. 



"QED" regime: Scalar potential disorder and 
flow toward weak coupling 



The "QED" regime occurs when the microscopic initial 
conditions are such that the (screened) scalar potential 
disorder fluctuations are much stronger than that of all 
other disorder types: g u ^> [see Eqs. (2.9)-(2.11), 
(2.26), and (2.14)]. Over some intermediate window of 
scaling, we would expect the RG to be dominated by g u 
and the Coulomb interaction strength w; to that end, we 
set all other disorder parameters to zero, and integrate 
the RG flow described by Eq. (4.12). 

As discussed in Sec. IV D 1, these flow equations pos- 
sess a repulsive fixed line, pictured as the dashed curve 
in Figs. 19 and 21. At the end of Sec. IVE, we noted 
that this line divides the g u -w plane into two portions; 
the portion joining the fixed line to the w-axis flows back 
to the clean, non-interacting Dirac fixed point, while its 
complement flows toward strong disorder and interaction 
coupling. (The portion that flows back to the clean Dirac 
fixed point is destabilized by the other disorder parame- 
ters for arbitrarily weak but non-zero disorder strength, 



24 



as shown in Appendix B.) The RG evolves toward weak 
(strong) coupling at longer length and lower energy scales 
if the "microscopic" disorder and interaction strengths g~ u 
and w satisfy g u < rjwf 3 {w) [g u > r)wf 3 (w)]. 

In the case of a flow toward weak coupling ("QED" 
regime), the Coulomb interaction w rapidly decreases 
toward zero; wc therefore take w <C 1 without loss of 
generality. One obtains the scaling behavior [using RG 
Eqs. (4.1a), (4.8), and (4.12), and Eqs. (5.3) and (5.4) 
for <r dc ] 



Ne 2 (rjwf 
2tt 2 lQg u 
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T J1 
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where 77 = 8/ttN and Tr 3> T is an arbitrary reference 
temperature. Setting N = 4, g u oc g u , and replacing 
w -> nr s /2 [Eq. (2.24)], we obtain Eq. (5.5) in Sec. V. 

The evolution towards weak coupling characteristic of 
the "QED" regime is represented via the schematic RG 
flow labeled (1) in Fig. 21. Although this flow approaches 
the clean, non-interacting Dirac fixed point [located by 
Eq. (4.16) in Sec. IV E], it is ultimately deflected to- 
ward the strong-coupling, "QCD" regime, labeled (3) in 
Fig. 21, for any non-zero, but arbitrarily small random 
vector potential (<7a) or random mass (g m ) disorder fluc- 
tuations. 



for T R > T. We recover Eq. (5.8) for N = 4. 

The evolution of the RG in this intermediate regime is 
represented by the curve labeled (2) in Fig. 21. 



3. "QCD" regime: Non-Abelian vector potential 
disorder and asymptotic flow toward strong coupling 

Finally, we consider the flow toward strong coupling 
for moderate to strong Coulomb interactions. This is 
the "QCD" regime dominated by the non-Abelian vector 
potential disorder parameters g& and <?A3 ■ As discussed 
in Sec. IV E, valley space SU(2) rotational symmetry is 
restored on average as the RG scales toward strong cou- 
pling; we may therefore take gA ~ gA3- 

Neglecting all disorder parameters except gA, we inte- 
grate Eqs. (4.11a), (4.11b), and (4. lid), and combine the 
results with Eqs. (4.8), (5.3), and (5.4). For w > 1, we 
obtain the following formula for the asymptotic scaling 
behavior of the conductivity: 
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for Tr > T. In this equation, 



2. Intermediate regime: Scalar potential disorder 
and flow toward strong coupling 

We next consider the case where again g u ^> Q^, 
but now g u > rjwf 3 {w), so that the RG flows toward 
strong coupling [Eq. (4.12)]. We further restrict w -C 1, 
since larger values of w induce a crossover to the "QCD" 
regime dominated by the vector potential disorder, de- 
scribed in the next subsection. 

With these assumptions, the intermediate scaling be- 
havior of the conductivity is given by 



a dc (T) 



Ne 2 
2^2 



9u \T 



(C2) 



1 

®9a 



1 - 



(C4) 



and W(z) is Lambert's function, i.e. solves the equation 
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For gA -C 1, Eq. (C3) reduces to 
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which is the same as Eq. (5.10) for V = 4 (77 = 2/w). 

The curves labeled (3) in Fig. 21 schematically indicate 
the flow toward strong coupling associated with the scal- 
ing prediction given by Eq. (C5) in this "QCD" regime; 
see also Figs. 15-18 in Sec. IV E. 
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